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Deep Ha Observations of NGC 253: a Very Extended and 
Possibly Declining Rotation Curve? 
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ABSTRACT 

This study presents a deep Ha kinematical analysis of the Sculptor Group galaxy 
NGC253. The Fabry-Perot data were taken with the 36-cm Marseille Telescope in La 
Silla, Chile, using an EMCCD detector. Typical emission measures of ~ 0.1 cm" 6 pc 
are reached. The observations allow the detection of the Diffuse Ionized Gas component 
through [Nil] emission at very large radii of 11.5', 12.8' and 19.0', on the receding side 
of the galaxy. No Ha emission is observed at radii larger than the neutral component 
(11.5'). The very extended rotation curve confirms previous results and shows signs of 
a significant decline, on the order of 30 per cent v max . Using the rotation data, mass 
models are constructed with and without the outer [Nil] data points, and similar 
results are found. The declining part of the rotation curve is very well modeled, and 
seems to be truly declining. 

Key words: galaxies: NGC 253 - galaxies: kinematics and dynamics - galaxies: ISM 
- instrumentation: interferometers - techniques: radial velocities. 



1 INTRODUCTION 

The Sculptor Group har bours an impressive amount of 
gas-rich late- type galaxies (|de Vaucouleursll 19591 ). including 
NGC 7793, NGC253, NGC 300, NGC 247 and NGC 55. 
Several other small gala xies are known to be gravitation- 
ally bound to the group l|Cote et al.lll997h . The galaxies are 
located at a distance estimated between 1.5 and 4 Mpc (Jer- 
jen, Freeman & Binggeli 1998; Karachentsev et al. 2003). Its 
proximity, as well as the fact that its galaxies are isolated, 
makes it an ideal ground for studying the kinematical prop- 
erties of disc galaxies. Moreover, the majority of its galaxies 
being late-type with no or almost no bulge, renders them 
excelle nt tools for studying the p roperties of dark matter 
haloes (|Carignan fc Freemanlll985l ). 

The literature contains many studies of the Sculptor 
Group galaxies. Among them, a deep Ha s tudy of NGC 
7793 was performed by iDicaire etafl 420081 1. This study, 
based on deep optical Fabry-Perot (FP) data, confirmed 
that NGC 7793 has a truly declining rotation curve. Deep 
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surveys of galaxies are essential to extract very extended 
rotation curves, and allow better constraints on the dark 
halo parameters. This study aims at studying the dynami- 
cal properties of another galaxy belonging to the Sculptor 
Group, NGC 253, also using deep Ha FP data. 

NGC 253 is a late-type spiral starburst galaxy. The cen- 
tral regions show significant activity, and have therefore been 
thoroughly studied. H53a and H72a observati ons showed 
the p resence of a counter-rotating central disc |Zhao et all 
l200ll ). Furthermore, observations of the H2 emission line 
suggest the presence of gas being ejected to distance s of 130 
parsecs above the disc jSugai. Davies fc WardllioO^ ). 

iPuche. Carignan fc van Gorkoml (|l99ll ) presented a 
kinematical study of the neutral gas component (H 1) using 
21-cm synthesis observations. A rotation curve extending to 
a maximum radius of 8.5 kpc (or 11.5') was derived, and 
showed that even at the largest radius, the rotation curve 
continues to rise with a maximum velocit y of 224 km s -1 . 
A long -slit Ha study has also been done bv lArnaboldi et all 
l|l995t} . 

The Diffuse Ionized Gas (DIG) is an important compo- 
nent of the interstellar medium that can extend to large 
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Table 1. Optical parameters of NGC 253 



Morphological type 
R.A (2000) 1 
Doc. (2000) 1 

Isophotal major diameter 1 , D25 (B) 
Holmbcrg radius 2 . i?HO (B) 
Exponential scale length 2 , a^ 1 (B) 
Axis ratio 2 (q = b/a) 
Inclination 2 (go = 0.15), i 
Total apparent B magnitude, B^ i5 
Total corrected apparent B magnitude, B^ 
Adopted distance (Mpc) 6 

Absolute B magnitude, Mg l 

Total blue luminosity (Mq — 5.43), Lb© 



CM3,4 



SAB(s)c 

00 h 47 m 33.1 s 

-25° 17' 18" 

27.5' 

17.5' 

2.4 kpc 

0.23 

79.8° 

8.05 

7.09 

2.58 

(l' = 0.75 kpc) 

-19.97 

14.4 x 10 9 



!RC3 

2 Puche et al. (19911 

internal absorption A(i) = 0.91 (RC3). 
4 Galactic extinction A B = 0.05 (RC3). 
5 Pence (1980 1 
6 Puche fc Carignan (1988 1 



radii. Studies have s hown that the gas is both warm, 
T ~ 6000 - 10000 K jRevnoldslll992l) and diffused, 



0.03 cm -3 (|Dopita fc Sutherland! 120051 '). However, the ori- 
gin of the DIG still remains controversial. One of the lead- 
ing theories implies ionizing photons escaping H n regions 
and reaching very large distan c es of several kiloparsecs. 
Irloopes. Walterbos fc Greenwaltl |l996) studied the proper- 
ties of this diffuse component in NGC 253. They measured 
the fraction of Ha luminosity from the DIG, and found it 
to be between 35 and 43 per cent. The distribution of the 
DIG also seemed to be concentrated near the Hn regions, 
extending away from them up to a maximum radius of 1 
kpc. 

The kine matical properties of the DIG in NGC 25 3 
were studied by Bland-H awthorn. Freeman fc Quinnl (|l997h . 
They found very extended Ha and [N n] emission on the 
south-west side of the galaxy at radii larger than that of the 
neutral component (Hi). The neutral component is detected 
out to 11.5', while the optical component extends to 15.0'. 
The last observational point of the optical rotation curve 
also suggests a decreasing rotation curve. Observing emis- 
sion at these radii could therefore have great implications, 
since it would allow the study of the distribution of the dark 
halo component at very large radii, and perhaps even allow 
the study of the entire distribution of matter. 

This paper presents deep optical FP data of the Sculp- 
tor Group galaxy NGC 253. A look for faint emission on both 
sides of the galaxy will be performed, as well as a search for 
emission beyond a radius of 1 5.0' in order to confirm the de - 
clining rotation curve seen bv lBland-Hawthorn et ail (119971 ). 
Section 2 will give a description of the FP observations and 
of the data reduction techniques used to extract the velocity 
field and the total integrated Ha emission map. It will also 
describe a method developed for the flux calibration of the 
Ha FP data. Section 3 shows the velocity field obtained in 
this study, and Section 4 presents the extracted kinemati- 
cal parameters and adopted rotation curve. The mass mod- 
els analysis is described in Section 5. The results are then 
analysed and discussed in Section 6 and Section 7 presents 
a summary and the conclusions. The optical parameters of 



NGC 253 are summarized in Tabl e [J A distance of 2.58 Mp c 
has been adopted for this study (|Puche fc Carignanlll9880 . 



2 FABRY-PEROT OBSERVATIONS AND 
DATA REDUCTION 

2.1 Observations 

The FP observations were obtained on the 36-cm Marseille 
Telescope at La Silla, Chile, from 2007 October 3 to 2007 
October 15. A commercial camera An dor iXon, equippe d 
with an EMCCD detector, was used (see lDaigle et alj|2004l '). 
All the observations were made with the same filter, camera 
and FP etalon. Dark frames were taken at the beginning 
and end of each night. Furthermore, data were taken to cal- 
ibrate the gain of the EMCCD. Flats with the etalon and 
Ha filter kept in place were also taken at twilight. The in- 
terference filter had a central wavelength of A c = 6561.0 A, 
a FWHM of 30 A, and a maximum transmission of ~ 80 
per cent at a temperature of 10° Celsius, see Fig. [4] The 
FP interference order was p = 765 at Ha, with a Free Spec- 
tral Range (FSR) of 8.58 A. The mean Finesse during the 
observations was ~16, with an etalon spectral resolution of 

~ 0.53 A (AA = ^ f^f^l ) ■ The data were scanned through- 
out the FSR over 40 channels of 0.21 A per channel. Wave- 
length calibration was performed using a neon lamp (Ne) 
at 6598.95 A. Several observations of the calibration cube 
were made during each night. For each observation of the 
calibration data, the etalon was scanned during one com- 
plete cycle (40 channels), with an exposure of 1 second per 
channel. The EMCCD has 512 x 512 pixels and a resolution 
of 2. 77" /pixel. Each field covers ~ 23.6' x 23.6'. 

The EMCCD has a quantum efficiency of ~ 90 per cent 
at Ha and the Electron Multiplying gain used yielded a 
read-out noise of a ~ 0.05 electron. In the Photon Counting 
Mode, the camera is operated at a higher frame rate and 
only one photon per pixel can be detected, which eliminates 
the excess noise factor. Nonetheless, in order to lower the 
effect of the Clock Induced Charges generated by the An- 
dor camera, the observations were done in Amplified Mode 
(15 s per image), where more than one photon per pixel 
per frame can be counted. This mode of operation causes 
the excess noise factor of the EMCCD to affect the pho- 
tometric quality of the data. More detail s can be found in 
iDaigle. Carignan fc Blais-Ouellettej (|2006t ) and lDaigle et al] 
(|2004l l2008i V 

NGC 253 was observed over three fields: centre, north- 
east and south-west. See Fig. [3] for outlines of the fields. 
Each field was exposed with an exposure time of 15 seconds 
per channel. The central field was exposed for 64 cycles, 
corresponding to 640 minutes or 10.7 hours. The north-east 
and south-west fields were each exposed for 38 cycles, cor- 
responding to 380 minutes or 6.3 hours. 



2.2 Data reduction 

2.2.1 Reduction steps 

Raw observational data consist of many data files taken for 
every cycle and channel. Each observational run consists of 
one to twelve cycles, and for each run, a wavelength-sorted 
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data cube can be extracted. The first step is to correct the 
data with the normalized flat image. In order to correct for 
the airmass dependance, which can vary considerably during 
an observing run, the correction must be applied prior to 
the creation of a wavelength-sorted data cube. These two 
corrections are applied at the beginning of the reduction 
procedure, for each and every channel of every cycle. The 
following steps were adopted for the reduction of the data: 

1. Correction of all data with the normalized flat image. 

2. Creation of a wavelength-sorted data cube corrected 
for airmass dependance, guiding shifts, cosmic rays, dark, 
and gain. 

3. Ghosts correctiorQ (see lEpinat et al.ll2008f ). 

4. Heliocentric velocity correction. 

5. Hanning sp ectral smoothing and sky subtraction (see 
iDaigle et afll200fj ). 

The data cubes were then summed for each field of NGC 
253, using the position of bright stars. The three fields were 
then mosaicked together, and corrected for the different ex- 
posure times using the computed flux of bright stars in each 
field. The final data cube has 1024 x 512 pixels. 

6. Spatial 5x5 binning procedure. 

7. Single gaussian fitting procedure. 

8. Correction for the filter dependance: division of the 
profile by the filter spectrum. 

9. Extraction of the velocity, velocity dispersion and Ha 
integrated maps. 

2.2.2 Lorentzian and gaussian profiles 

The intensity weighted method, which is often used to ex- 
tract the velocity of galaxies through spectral-line profiles, 
calculates the barycentre of the spe ctral line, and is appro- 
priate for cases of low S/N data (see lAmramlll99ll . for more 
details). However, the shape of the emission lines observed 
should be the result of a convolution between the device 
function, an Airy function, and the emitted spectral line, 
which is assumed to be gaussian. Thermal agitation of the 
gas dominates, and is responsible for the gaussian shape of 
the profile. The observed profile has therefore been approx- 
imated by a Lorentzian or a gaussian profile. The gaussian 
profile yielded a mean reduced x 2 r °f 1-37 for the entire 
galaxy, while the Lorentz profile yielded a mean x 2 r °f 2.67. 
The gaussian fitting procedure was therefore adopted. 

2.3 Ha flux calibration 

The Ha flux calibration of the data was done during the 
reduction process. In order to carry out the calibration, an 
extinction coefficient as well as a conversion factor (FAC) 
must be determined as defined by Eq. [T] and Eq. [2] where 

1 Ghosts are reflections of bright Hn regions in the field. Spa- 
tially, the ghosts are located symmetrically opposite with respect 
to the optical axis and they will be at the same wavelength as the 
primary image. They can thus be well identified, both in position 
and in velocit y. A ghost co r rectio n was applied using the method 
developed by lEpinat et al.l II2008T ) . Some residuals of ghosts were 
seen, but these fell far away from the galaxy. They also did not 
fall in the pie-shape regions used to study the DIG (sec Section 
3 and 4). 



Flux is the corresponding flux value, X is the airmass at the 
time of exposure, Cobs is the observed number of counts, 
and Ctmc is the true number of counts. In order to alleviate 
the total number of parameters to adjust, the zero point 
parameter, B, and the initial conversion factor, / ac , were 
combined together to obtain the conversion factor (FAC). 

FluX = / ac • Ctruc 

= /ac • W ^ X + B ^ 5 ■ C ohs 

= 10 (A ' X)/2 - 5 • C obs ■ /ac • 10 (S/2 - 5) (1) 

Flux = i ( A ' x >/ 2 ' 5 . C obs . (FAC) (2) 

For this study, the value of the extinction coefficient pub- 
lished by ESO was used. The val ue for Ha at the La Silla 
site is 0.07 mag/airmass (see also lBurki et al.lll995T ). 

The next step was to extract the conversion factor 
between the number of counts corrected for airmass and 
the corresponding flux. This was accomplished once all the 
reduction steps were performed, and the total number of 
counts was obtained, since it allows a better constraint 
on the factor determined. The continuum subtracted Ha 
map from the SHASSA (Southern Ha Sky Survey) cata- 
logue was used to flux calibrate the data (|Gaustad et al.l 
l200ll . |http: / /amundsen.swarthmore.edu/l ) . This catalogue 
has units of deci-Rayleigh, dR (1R = lO* 3 photons cm~ 2 s _1 
sr _1 or 0.6 x 10 -17 ergs cm~ 2 s _1 arcsec -2 or to an Emis- 
sion Measure, EM, of ~ 2 cm -6 pc). Moreover, the SHASSA 
image was corrected for [N n] (AA 6548 A, 6583 A) contam- 
ination. For NGC 253, the intensity of the contamination 
was considered to be 0.48 times that of the Ha intensity 
l|Kennicutt et al.ll2008T I. Afterwards, a linear regression was 
applied between the counts observed from this study and 
those of the SHASSA catalogue for individual bright Hn re- 
gions in NGC 253. Fig. [1] compares the calibrated data from 
this study with those of the SHASSA image for bright Hn 
regions. 



3 VELOCITY FIELD 

The velocity field was extracted using three distinct meth- 
ods. The first followed the reduction steps outlined in Sec- 
tion 2.3.1, and the resulting velocity field is presented in 
the top-left image of Fig. [5] The corresponding velocity 
dispersion map and Ha integrated map are also shown. 
As expected from an inspection of the velocity dispersion 
map, the velocity field contains perturbations, notably in 
the inter-arm region of the northern side. In order to deter- 
mine the global kinematical parameters, a second method 
using a gaussian smoothing with w\ = 55" (20 pixels) was 
used to extract the velocity field, where W\ is the width 
in " of the smoothing function (see bottom-right image of 
Fig. [2J. In order to study the very diffuse emission at ex- 
tended radii, a third method was used to extract the veloc- 
ity field. A special binn i ng pro cedure similar to that used by 
iBland-Hawthorn et al.l (|l997h was applied. The smoothing 
function has the shape of an isoscele triangle (pie-shape), 
centred around the major axis of the galaxy and about its 
dynamical center. For a given opening angle, different re- 
gions corresponding to different radii intervals were binned 
together. An opening angle of ~ 15° centred around the ma- 
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Figure 2. Kinomatical maps of NGC 253 obtained after a 5 pixels X 5 pixels square binning and a single Gaussian fitting procedure, 
unless otherwise specified. Top-left: Ha velocity field. Top-right: Dispersion map. Bottom-left: Flux calibrated integrated map obtained 
after correction for [Nil] contamination (scaled in 1/6 Raylcigh or 10~ 18 ergs cm -2 s _1 arcsec -2 ). Bottom-right: Ha velocity field with 
a gaussian spatial smoothing of w\ = 55". 




5 10 15 20 

Calculated Ha Flux (10~ ,s erg cm"' s" 1 arcsec"*) 

Figure 1. Comparison be tween the fluxe s obta ined for NGC 253 
in this study, and those of lGaustad et all J200lf) . Each point rep- 
resents a bright Hn region in NGC 253. The conversion factor 
was determined to be 4576 (1 count = 4576 dR, where 1/6R or 
an EM of 2/6 cm -6 pc is equivalent to 10 — 18 ergs cm~ 2 s _1 
arcsec -2 ). The continuous line has a slope set to unity. 



jor axis was chosen, with the binning intervals growing with 
radius (see Fig. [3]). 

The central regions of NGC 253 are very bright. It is 
therefore possible that this bright source causes the appear- 
ance of parasitic emission in the field, similar to ghosts. The 
light from the nucleus is a strong continuum source. When 
passing through the filter, diffusion of the light can be ob- 
served. These ghosts will appear as rings in the velocity map, 
on top of the velocity field of the galaxy. The following sec- 
tion describes in detail this parasitic emission, as well as the 
method used to distinguish it from the diffuse emission of 
the galaxy. 



4 ISOLATING THE EXTENDED DIFFUSE 
EMISSION 

4.1 Diffuse Ionised Gas 

The term DIG refers to a component of the i nterstellar 
medi um that is both warm, T ~ 6000 - 10000 K dRevnolds 
19921) , and diffused, n B ~ 0.03 cm" 3 (|Dopita fc Sutherland 
2003). The DIG has low values of EM, since the emission 



measure depends on the square of the density. This com- 
ponent can be studied by separating its emission from the 
bright, high EM, ionised regions in a galaxy, such as Hll 
regions. Some have simply separated this component by ap- 
plying a cut off: any emission below a certain EM, such as 
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Figure 3. Top-left: Velocity field of NGC 253 obtained after a 5 X 5 square binning procedure, a pie-shaped binning procedure, and 
a single gaussian adjustment to the emission profiles. The last three binned intervals (black contours) are the result of [Nil] emission. 
Outlines of the three fields arc also plotted in large squares: black for centre field, blue for north-east field and red for south-west field. 
Top-right: Close up of same velocity field, but without any masking applied - the rings of the parasitic emission embedded in the velocity 
field are clearly seen. The black annuli contours show the regions binned to extract the typical parasitic spectra. Bottom: Spectra of 
the last three binned intervals of the pic-shape in order of radius (from left to right). The star-like points are for the spectra of the 
DIG and the continuous black profile are for the corresponding typical parasitic spectra. The blue vertical line represents the location 
of the Ho line, the red vertical lines represent those of the [N II] emission, and the grey regions give an indication of the location of the 
contaminated areas. 



~ 80 pc cm 6 (value taken from iFerguson et al.lll996l ). was 
considered to be from the DIG. 

However, iHoopes et al.l (|l996l ) argue that this simple 
cut off method is not adequate for highly inclined galaxies 
such as NGC 253, and that it can result in underestimating 
the total DIG component. In their study, they developed 
a special median filtered masking technique to isolate the 
DIG e mission from their NGC 253 images, see lHoopes et al] 
|l996t ) for more details. 

Nonetheless, the goal of this study is not to separate the 
DIG emission from the bright ionised regions. The goal is to 
obtain a very deep image of the galaxy in order to study 
the kinematics at very extended radii. This does not require 
isolating the DIG component, since any emission whether 
with high or low EM can be used to obtain kinematical 
measures. 

For the extended ionised emission detected beyond ~ 
11', we use the term DIG. It is expected that at large radii, 
the EM of the gas drops significantly. For radii beyond ~ 
11' in NGC 253, the EMs are well below a few pc cm -6 , 
which qualifies it has DIG regardless of the inclination of 



the galaxy. In the following section, we explain how this 
very low EM gas was detected, and how its emission was 
isolated from the parasitic emission mentioned in Section 3. 



4.2 Isolating the extended DIG emission from the 
parasitic emission 

The parasitic emission seen in the data of this paper is 
thought to originate from diffusion of the bright nucleus 
light when passing through the filter. The same parasitic 
emission was seen on other observing runs, when the targets 
were bright point-like planetary nebulae. Since the nucleus of 
NGC 253 is particularly bright and point-like, the parasitic 
emission seen could therefore have the same cause which is 
thought to be diffusion of light through the filter. 

The top-right panel of Fig. [3] shows a close up of the 
raw velocity field (no masking applied), where the parasitic 
emission is clearly seen as concentric circular rings embedded 
in the velocity field of NGC 253. The centre of the parasitic 
rings coincides with the FP optical centre, and so rings were 
seen for each field (centre, north-east, south-west). 



© 2009 RAS, MNRAS 000.ITlfl5l 



6 J. Hlavacek-Larrondo et al. 



In Fig. [2] the velocity field of the galaxy extends up to a 
radius of ~ 11'. The intensity of the emission lines that gave 
forth to this velocity field are on average 6 times stronger 
than the intensity of the parasitic emission. Thus, for regions 
of strong Ha emission, i.e. up to 11', the parasitic emission 
is not problematic since the emission lines of the gas can be 
easily identified. The spectrum of each pixel in this velocity 
field was also individually examined, and all data that did 
not clearly show a strong emission line were discarded. 

For very low EMs, i.e. beyond a radius of 11', the line 
intensity can be as low as the intensity of the parasitic emis- 
sion. Diffuse emission can therefore be hidden in the signal 
of the parasitic emission. However, Fig. shows that the 
parasitic emission is not random noise; it has a well de- 
fined spectra which gives circular rings in the velocity field. 
Therefore, a typical spectrum of the parasitic emission can 
be extracted for each radius interval, and used to identify if 
the emission seen is from the extended DIG or the parasitic 
emission. 

The method used to isolate the extended DIG emis- 
sion is as follows, and is similar as the one used by Bland- 
Hawthorn et al. (1997) to identify the DIG emission in their 
study. This method is applied to all radii beyond 11', where 
diffuse emission is seen. The DIG emission was extracted 
using a pie-shape binning procedure centred around the dy- 
namical centre of the galaxy (see pie-shape in top-left of Fig. 
[3]). The spectra of the parasitic emission was extracted us- 
ing another pie-shape binning procedure, this time centred 
around the centre of parasitic rings (see pie-shape outlined 
by black annuli in top-right of Fig. [3]). The opening angle 
was chosen to be large enough so that a typical signal for 
the parasitic emission could be extracted, while smoothing 
out any emission from the galaxjjf). The radii intervals were 
chosen to be the same as the first binning procedure. Any 
emission profiles detected in the first procedure (pie-shape 
centred on galaxy) and not in the second (pie-shape centred 
on parasitic rings) should be the result of emission from the 
DIG. 

The bins with more than a 3<r detection are illustrated 
in the bottom panels of Fig. [3] These all lie on the receding 
side of the galaxy. They are the last three intervals shown 
with black contours in the pie-shaped velocity field and are 
located at 11.5'±0.6', 12.8'±0.7', 19.0'±0.9' (top-left of Fig. 
[3| . In these three bins, the bottom panels of Fig.0show that 
the DIG spectra (star-like points) have an emission line (blue 
or red line) not seen in the corresponding typical parasitic 
spectra (black curve). 

Different opening angles were tested for the pie-shape 
centred on the parasitic rings, but they always gave the 
same results. First, that the extracted spectra of the par- 
asitic rings did not vary significantly, which means that the 
parasitic emission is not just random noise and that a typ- 
ical spectra can be extracted. Second, when compared to 
the parasitic spectra, the DIG spectra always showed a 3cr 
detection for the same four emission lines mentioned in the 
previous paragraph. These emission lines lie along the major 



2 This gave an opening angle of about 140° , except for the first 
radii bin (r = 11.5' ± 0.6') where a much smaller opening angle 
had to be chosen. An opening angle of 140° would have overlapped 
significantly on emission from the galaxy. 



axis of the galaxy, are not from the parasitic emission, and 
therefore can be interpreted as emission from the DIG. 

The integrated flux of the first, blue, line (r = 11.5' ± 
0.6') gives an EM of ~ 0.14 cm" 6 pc, and those with the red 
line (r = 11.5' ± 0.6', 12.8' ± 0.7', 19.0' ± 0.9') give an EM 
of ~ 0.10 cm" 6 pc. Table [3] shows these results. The EMs 
were calculated by subtracting the contamination spectra 
from the DIG profiles, then correcting for the filter depen- 
dance and integrating the profiles. The la upper limit for 
the EM is ~ 0.03 cm" 6 pc. This was calculated by taking 
the standard deviation (a) between the spectrum resulting 
from the pie-shape centred on the galaxy, and that centred 
on the parasitic emission, for all the radii bins beyond 11'. 
For the three radii bins where DIG was detected, the emis- 
sion line profile was omitted before calculating a. From the 
averaged ct, and a typical width for the extended DIG emis- 
sion lines, the la upper limit for the EMs was found to be 
~ 0.03 cm -6 pc. Hence, the three emission lines seen have 
EMs on the order of 3a. Any emission below this limit was 
not considered as DIG emission. The error of the emission 
line fluxes is thought to be on the order of a, therefore of 
~ 30 per cent. 

4.3 Nature of the emission lines 

The extended DIG emission line found in the first radii in- 
terval (r ~ 11.5') and identified with a blue line in Fig 3 
should be the Ha line at A ~ 6572 A, according to a ve- 
locity pattern continuity argument. Due to the periodicity 
of FP spectrum, in this figure, the wings of this line are 
seen on both extremity of the spectrum. This is the phase 
jump effect. Phase jumps are common in FPs when the filter 
bandpass is larger than the FSR of the FP. 

More precisely, a phase jump, also known as 
interordre confusion of neighbouring FP orders (see 
iBland-Hawthorn et ailll992l ; lDaigle et al.ll2006r i. appears be- 
cause the interference pattern of the FP will transmit (for 
a given spacing of the plates) all wavelengths A such that 
A = Ao ±n x FSR where n is an integer value, providing that 
these wavelengths are also transmited by the intereference 
filter used for selecting the wavelength range to be scanned 
with the FP and to limit the number of phase jumps (see 
Fig- HI- Essentially, the FP acts like a secondary filter, with 
a transmission peak occurring every 8.6 A. In other words, 
light first travels through the interference filter, which af- 
fects the intensity of the light, depending on the transmis- 
sion probability (Fig. [4}. Then, light travels through the FP, 
where phase jumps can occur, but do not affect the inten- 
sity of the light. Therefore, the probability of having phase 
jumps of a specific emission line depends only on whether 
the probability of it passing through the interference filter 
is high enough. 

In this study, we have choosen the most adequate in- 
terference filter to cover the wavelengh range of the Ha line 
in NGC 253, given its systemic velocity and we study the 
FP interference order corresponding to the wavelength range 
between ~ 6564 A and ~ 6572 A, covering Ha velocities 
between 40 and 450 km s _1 (see Section 2.1). Hence, the 
spectrum is plotted between 6564 A and 6572 A, but there 
could be emission from neighbooring free spectral ranges: 
for instance, an emission line detected at A = 6570 A could 
be emitted at A = 6570 ± n x 8.6 A. 
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Figure 4. Interference filter transmission profile (black curve), 
shown in per cent. The positions of the redshifted (red) and the 
blueshifted (blue) [Nll]AA 6548, 6583A and Ha lines arc also 
shown. NGC 253 has a systemic radial velocity of ~ 241 km s _1 . 
The redshifted velocities were calculated using a radial velocity 
of ~ 450 km s _1 and the blueshifted velocities were calculated 
using a radial velocity of ~ 40 km s" 1 . As can be seen, the red- 
shifted [Nil] A 6548A is almost perfectly transmitted, much more 
than the redshifted Ha line. The FSR of the FP (8.6 A) is also 
shown in the upper-left corner. 



The DIG emission lines identified by a red line in Fig. 
[3] at A ~ 6566.4 A cannot be from Ha since the velocities 
implied are not consistent with the velocity pattern of the 
galaxy. They are rather emission from the [Nil] line (A 6548 
A) redshifted, and are probably observed around this wave- 
length due to a phase jump. 

Fig. 0] shows the filter transmission profile, and the lo- 
cations of the blueshifted and redshifted emission lines of 
[Nn]AA 6548, 6583A and Ha lines. This figure clearly shows 
that on the receding side of the galaxy (redshifted), the 
[Nn] A 6548 A line is more likely to be transmitted than the 
Ha line. The line could be emitted at 6566.4 A or at 6557.8 
A (in case of a phase jump) or even at 6549.2 A (in case of 
multiple phase jumps). 

Due to the interference filter transmission, this line is 
most likely emitted at 6566.4 A ( 75 per cent) or at 6557.8 
A (78 per cent). In the first case, the line we see at 6566.4 
A is an extremely redshifted [N n] line with a radial veloc- 
ity of ~ 850 km s _1 . This radial velocity is not consistent 
with the galaxy's velocity field, but could be interpreted as 
an accretion phenomenon. We examine this possibility in 
Section 7.2.1. In the second case, the velocity (~ 400 km 
s _1 ) is compatible with a declining rotation curve. We use 
this assumption in Sections 5 and 6. Why no extended DIG 
emission was seen on the approaching side of the galaxy, as 
well as why no Ha was seen on the receding side beyond 
11.5' will be looked at in Section 7.2.1. 



5 KINEMATIC AL PARAMETERS AND 
ROTATION CURVE 

In order to derive the rotation curve of NGC 253 from the 
analysis of the velocity field, the GIPSY task ROTCUR 
(|Begeman|[l989h was used. The objective is to extract the 
orientation parameters (dynamical centre (xo, yo), systemic 
velocity V sy s, inclination i and position angle PA) that best 
reproduce the observed velocity field. As mentioned ear- 
lier, the velocity field smoothed with a gaussian function 
of w\ = 55" is used. The parameters are found by a least 
square method between the observed velocity field and the 
tilted rings derived by ROTCUR. 

There are three main steps with the GIPSY program. 
Since they are correlated, the dynamical centre and systemic 
velocity are found first by keeping i and PA fixed at their 
initial value s. The initial value s are the ones derived by the 
H I study of lPuche" et all (|l99lh . Data in an opening angle of 
65° around the minor axis are excluded in order to minimize 
errors due to deprojection effects. A cosine-square weighting 
function is then used on the rest of the data, maximizing the 
weight of the points around the major axis. From this, the 
dynamical centre is found to be located at a = 00 : 47 : 35.2 
and 8 = —25 : 17 : 07.6, in the World Coordinate System. 
A systemic velocity of 241.2 ± 5.4 km s _1 is also found. 
Th is value is compati ble with the values derived from H I 



by IPuche et all (Il99lf) (~ 2 45 km s" 1 ) and from Ha by 
iBland-Hawthorn et all (|l997f ) (243 ± 9 km s" 1 ). 

The second step seeks to find the inclination i and the 
position angle PA, while keeping the dynamical centre and 
systemic velocity fixed at their extracted values. The initial 
values for the inclination i and p osition angle PA ar e again 
those derived by the Hi study of lPuche et al.1 (|l99ll ). Sepa- 
rate solutions are obtained for the receding and approaching 
side. A mean value for the inclination i of 76.4 ± 3.9° and 
for the position angle PA of 230.1 ± 2.1° are derived be- 
tween a radius of V and 10'. ROTCUR could not converge 
for a small number of radii, and the value of the parameters 
at these radii were extrapolated with the values of nearby 
data. The gaussian smoothed velocity field extends to a ra- 
dius of 10', while the original velocity field extends to 11.5', 
although the later is very clumpy beyond 10'. In order to de- 
rive the rotation curve from the original velocity field using 
the derived kinematical parameters, extrapolated values of 
the parameters were added for radii between 10' and 11.5'. 
The top figures of Fig. [5] show the adopted values for the in- 
clination i and position angle PA. The results a re consistent 
with t he parameters derived by the H I study of lPuche et al.1 
(|l99ll ), and also suggest the presence of a warp in the outer 
regions as seen through the decrease of the position angle 
at r ~ 5'. In the inner 5' region, a variation of the position 
angle is also seen. However, this can be attributed to the 
perturbations of the nucl ear region since NG C 253 is known 
to be a starburst galaxy. IPuche et al.1 |l991) find f° r the H I 
data a mean inclination i of ~ 72° and a mean position angle 
PA of ~ 229°. From the adopted values of the inclination 
i and position angle PA varying with radius, the rotation 
curve of NGC 253 could be extracted from the original ve- 
locity field ( without the gaussia n smoothing) . The warp seen 
in NGC 253 (|Puche et al.lll99ll ) can be accounted for by let- 
ting the inclination and the position angle vary with radius. 
The bottom-left image of Fig.|6]illustrates the rotation curve 
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Figure 5. Velocity field of NGC 253, in which the pic shape used 
to identify the diffuse emission is shown. Contours of the near- 
UV GALEX image in a logarithm scale are shown. The outermost 
contours was taken as that with 0.005 counts/second/pixcl. 



derived, along with the solutions found for the approaching 
and receding side. Finally, the bottom-right image gives the 
adopted rotation curve of the galaxy, with the proper error 
bars. The error bars correspond to the biggest difference be- 
tween the rotation curve of the whole galaxy and that of one 
of the sides, or, if larger, the intrinsic error determined by 
the ROTCUR task. 

The velocities determined for the DIG component from 
the previous section are shown in the bottom-right panel 
of Fig. [5] (diamonds) . The kinematical parameters used to 
derive the rotation velocity from the line of sight velocity 
were those at the maximum radius of 690". The error bars 
in the abscissae direction were chosen as the half width of 
the radii interval from which the emission line was extracted, 
while the one in the vertical direction was carefully chosen 
as the one resulting from the dispersion of the emission line 
since the largest uncertainty is on the precise location of 
the profile, rather than on the variation in the kinematical 
parameters. 

In order to derive the velocities of the DIG, we assume 
that circular motions remain at large radii. Non-circular mo- 
tions have been studied in the central regions of galaxies, 
where phenomena such as bars have been found, but it re- 
mains difficult to study these motions in the outskirts of 
galaxies. Dee p surveys of matter at large radii are needed. 
The study bv lDicaire et al.l (|200ct ) found an extended Ha ro- 
tation curve for NGC 7793, another galaxy of the Sculptor 
group. The Ha emission extended out to the H I, and showed 
a rotation curve consistent with that of the H I. When deriv- 
ing the kinematical parameters of NGC 253, no significant 
non-circular motions were seen at r > 7'. For NGC 253, we 
examined the surface-brightness isophotes of the near and 
far ultra-violet (UV) maps, taken from the Galaxy Evolu- 
tion Explorer (GALEX) survey. On the receding side of the 
galaxy, the isophotes extended out to a radius of about 13'. 
This wavelength traces hot stars, which could be a source 
of ionization and could explain our extended observations. 



Table 2. Adopted Ha rotation curve 



radvus 


T/ 

V rot 


A 1/ 
^ V rot 


I'll lit us 


1/ 

V rot 


£i V rot 


(axcsec) 


(km s~ 


(km s~ 1 ) 


(arcsGc) 


(km s~ 


(kin s~ 


60 


123.96 


35.66 


390 


183.09 


19.21 


90 


124.54 


9.97 


120 


182.45 


22.31 


120 


164.25 


4.97 


450 


185.27 


20.80 


150 


189.73 


6.71 


ISO 


198.82 


9.74 


180 


201.20 


3.50 


510 


199.55 


4.01 


210 


197.36 


2.82 


540 


203.52 


1.54 


240 


192.80 


10.79 


570 


202.81 


10.88 


270 


187.55 


13.65 


000 


212.54 


9.65 


300 


170.12 


30.85 


630 


212.72 


11.67 


330 


160.46 


36.60 


BOO 


208.00 


13.81 


360 


168.38 


32.31 


090 


197.72 


18.30 



The isophotes showed distortions in the inner regions, but 
were symmetric out to 13' (see the near-UV contours shown 
in Fig. [5|. This suggest, that at least for the two inner bins 
where we detected diffuse gas (r = 11.5', 12.8'), we can 
treat the rotation as been circular. For the outermost point 
(r — 19.0'), our observations do not allow us to make a 
definitive answer. 

A deep op ti cal im age of NGC 253 was taken by 
Malin fc Hadlevl (119971 ). and is presented in Fig. 2 of 
Boomsma et al.l (|2005r ). This image supports the idea that 
the stellar disc of NGC 253 extends well beyond the H I 
disc, but the image also suggests that there are some distur- 
bances on the south-west side of the galaxy. This could con- 
tradict the results regarding the Galex isophotes. However, 
w e prefer to keep with the assu mption of circular orbits as 
in lBland-Hawthorn et all (| 19971 ) 's paper, and see what kind 
of results we obtain. 

The adopted values of the rotation velocity and the 
associated errors are given in Table [5] and Table [3] In or- 
der to compare our resul ts, with the previou s ones, the H I 
rotation curve found by IPuche et alj (|l99ll ) is illustrated 
in Fig. H3 as well as t he ext ended data points found by 
iBland-Hawthorn et alj (|l997l ). where the inclination and 
systemic velocity of this study were used. For the later, 
see Table 3. Finally, the Ha rotation curve derived by 
lArnaboldi et~aH (|l995l ). using long-slit observations, is also 
illustrated in the figure. The latter study did not correct the 
rotation curve for the inclination or the PA, and therefore 
did not consider the variation in the kinematical parameters 
with radius. Deviations are therefore expected. Their curve, 
plotted in Fig. HJ] has been corrected for inclination using 
the mean value derived for our study. 



6 MASS MODELS 

In order to determine the impact of the DIG component ve- 
locity points at large radii, two mass models analysis have 
been obtained for NGC 253. The first has been performed 
on the data taken without the three points of the DIG 
while the second included those extended velocities. Fur- 
thermore, in order to compare our resu lts with those found 
for the Hi data bv lPuche et alj ( |l99lf ) . the analysis of the 
mass distribution must be the same. Hence, another mass 
model analy s is has been obtained using the data taken from 
IPuche et alj (|l99ll ). A dark halo represented by an isother- 
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Figure 6. Top: Kinematical parameters of NGC 253 derived from the gaussian smoothed (u>\ = 55") velocity field for both sides 
(circles) and separately for the receding (red triangles) and for the approaching (blue triangles) sides. The non-filled symbols represent 
the extrapolated parameters, see Section 4 for more details. Bottom-left: Rotation curve derived using the kinematical parameters of 
the top figures as a function of radius, with a systemic velocity of Vaya = 241.2 km s _1 , and the non gaussian smoothed velocity field. 
The rotation curve obtained using both sides at once is represented by the continuous line, along with the derived error bars, while the 
receding (red) and approaching (blue) sides are represented by the triangles. Bottom-right: Final ad opted rotation curv e, with the proper 



19911). and the double- 



error bars as explained in the text (continuo us line). The slashed lin e represents the H I data from iPuche et al.l I 

dot sl ashed line represents the Ho data from I Arn aboldi ct al. ( 1995]). The colored points illustrate the results of ] 

(1997) for [Nil] (orange squares) and for Ha (green circles) beyond the Hi data, as well as those derived from this study (blue and red 
diamonds). 
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Table 3. Derived DIG rotation velocities and those of Bland- 
Hawthorn et al. (1997) 

This Study 



radius 


Ar 


Vrot 


AVrot 


EM 


(arcsec) 


(arcsec) 


(km s" 1 ) 


(km s" 1 ) 


cm -6 pc 


690" 


36 


200.51 
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0.14 
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42 
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51 
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630 6 


30 
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8.3 


0.18 


720" 


25 


223.2 


8.0 


0.16 


810 6 


20 


229.0 


8.3 


0.09 


900" 


15 


200.8 


25.0 


0.08 



a Ha 
b [Nil] 



mal sphere is used. The density of the isothermal sphere 
has a fiat slope in the central regions. Observationally, the 
data favour models with flat slopes rather than those with 
cuspy c ores such as the Navarro, Frenk fc White (NFW) pro- 
file fe.g.lBlais-Ouellette et alj|200ll:lde Blok fc Bosmall2002l : 



Ide Blok et al.ll2003l ; iKassin et alj|2006j ) 



Three parameters must be specified: two for the dark 
halo, namely a one-dimentional velocity dispersion a and 
a core radius r c , related together by the central density 
po = 9a 2 /4irGr c 2 , and the mass-to- light ratio (M /L)* o f the 
stellar disc. The papers bv lCarignan fc Freeman! (1 19851 ) and 
ICarignanl (| 19851 ). show in greater detail the method used. A 
least square fit is applied to determine the three parameters. 
In order to compare with the H I study, the stellar disc is de- 
sc ribed by the surfac e brightness profile in the B band, as 
in IPuche et all (|l99ll ). The gaseous component is also taken 
as the one derived in the H I study. The latter includes the 
4/3 multiplication factor to account for the H e contribution. 

The mass model is designed to put higher weight on 
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Figure 7. Mass model analysis for the adopted Hi rotation 
curve (star-l i ke po ints) of NGC 253. The data were taken from 
iPuche et al.l d 199lt) . An isothermal sphere was used for the dark 
halo. A first analysis was performed using a maximum stellar disc 
(top), and a second using a sub- maximum stellar disc (bottom). 
The stellar disc is plotted with a small-dashed line, the gaseous 
component with a dot-dashed line, the halo with a long-dashed 
line and the total adjusted rotation curve with the continuous 
line. The parameters extracted are in Table. l4l 



the data points that have small errors. In the case of the 
Hi, the data at intermediate radii, between 3 and 5.5 kpc, 
have much smaller errors than the central and outer points. 
On the other hand, in the case of the Ha data, the mass 
model tended to fit as best as possible the central areas 
which had the smallest errors. This yielded very different 
results. Hence, in order to better compare the Hi and the 
Ha data, equal weight was given to all data points. 

First, a maximum disc model with no halo component 
was obtained (see top of Fig. [Jj Fig. [8j Fig. [9] and Ta- 
ble 2J). Second, a sub-maximum disc was obtained. Here, 
(M/L)* is fixed to the value that allows the maximum ve- 
locity of the stellar disc to be 0.6 times the maximum ve- 
locity of the rotation curve (Hi or Ha) as suggested by 
the study of ICourteau fc Rixl (119991 ). Indeed, the (M/L)* 
determined from the maximum disc is usually too high to 
agree with the predict ions of the stellar populations models 
|Bell fc de Jondl200ll ). See bottom of Fig. [JJ Fig. [8j Fig. [9] 
and Tabled 
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Figure 8. Mass model analysis for the adopted Ha rotation 
curve from this study (star-like points) of NGC 253, without the 
DIG [Nil] points. Same caption as Fig. [7J with a maximum stel- 
lar disc (t op) and a sub-max imum stellar disc (bottom). The Hi 
data from lPuche et aL 1 < 19911) are also plotted with the diamond- 
shaped points. 



Table 4. Mass models for NGC 253 





Hi RC£ 


36-cm Ha RC 
without DIG 


36-cm Ha RC 
with DIG 


Max Disc 








(M/L B )*Ji 


12.0 


12.9 


12.5 


Dark halo^ 


none 


none 


none 


Sub-Max Disc 








(M/£b)* 


5.0 


5.0 


5.0 


Dark halo 


r c = 2.3 
cr = 2.72 
po = 0.286 


7-c = 1.0 
a = 3.04 
po = 1-579 


r c = 1.0 
a = 2.99 
po = 1.564 



a Puche et al. (1991) 
^ in units of Mq /Lq 

c r c in kpc, <r in km s -1 and po in Mq pc — 3 
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Figure 9. Mass model analysis for the adopted Ha rotation 
curve with the DIG [N o] + Ha points (star-like points) of NGC 
253. Same caption as Fig. [7] with a maximum stellar disc (top) 
and a sub-maximu m stellar disc (bottom). The Hi data from 
iPuche et al. (1991) are also plotted with the diamond-shaped 
points. 



7 ANALYSIS OF THE RESULTS 
7.1 Nuclear region 

NGC 253 is a starburst spiral galaxy. The central regions 
have been thoro ughly studied. S ome suggest that NGC 253 
is also an AGN (|Ptak et aL 19971) and has extensive galactic 
wind s (see a review by IVeilleux. Cecil fc Bland-Hawthorn! 



120051 ). The paper by Heesen et al. (2009) also finds that the 
poloidal component of the magnetic field sh ows a prominent 
X-sha ped structure centred on the nucleus. iBoomsma et al.l 
|2005l ) also find extra-planar H 1 on the north-east half of the 
galaxy, in the shape of a half-ring structure and plumes, with 
a smaller rotation velocity than the disc. It is suggested that 
the origin of the extra-planar H I is related to the starburst 
activity and the active star formation in the disc. 

The present study finds high Ha emission and high 
velocity dispersions in the nuclear region, which are ex- 
pected from the strong nuclear activity. Higher-resolution 
dat a have also seen a large dispersion in the velocity field 
(see lZhao et al.|[200l[ ). which suggest that the one found in 
this study is not due to beam smearing. A small bar-like 
structure is also seen in the velocity field, lagging in velocity 
by at least 50 km s _1 (see top-left image of Fig. [2] the bar- 
like structure in the nuclear region appears more blue-shifted 



than its surroundings). The structure also seems slightly in- 
clined with respect to the major axis. The dimensions are 
approximately 24" by 88", i.e. 0.3 kpc by 1.1 kpc. The strong 
activity in the nuclear region of the galaxy could account 
for thes e observations, inclu ding the presence of a possi- 
ble bar (jPaglione et al.ll2004l) o r a po ssible counter-rotating 
circumnuclear disc ( Zhao et alj|200lh . An inclined counter- 
rotating disc could be overlaying the bulk circular rotat- 
ing disc a nd account for the lagging velocit ies a nd strong 
disper sion. iRubin. Graham fc Kem icv (1992j) and lRix et al.l 
|l992l) have shown for the early-type galaxy NGC 4550, that 
a natural explanation for a decrease in velocities and increase 
in dispersion could be the presence of two counter-rotating 
stellar components, as indicated by the observed double- 
peaked line profiles. Although NGC 253 is not an early-type 
galaxy, but a late-type, this explanation could also be ap- 
propriate for its gaseous component. 

In Fig. [6] the Ha rotation curve derived from this study 
is seen to ag r ee we ll with the Hi rotation curve derived by 
IPuche et al.l (|l99ll ). except for the inner regions. Here, the 
Ha slope is much steeper. This can be understood since the 
Hi data is of much lower resolution (~ 68") compared to 
the Ha data (~ 2.8"), and beam smearing is expected which 
causes the inner slope to be less steep. Hence, an important 
advantage is that Ha data probes more accurately the inner 
regions of the rotation curve. 



7.2 Diffuse Ionized Gas 

7.2.1 [Nil] and Ha emission 

This study succeeded in detecting the [N 11] (A 6548 A) emis- 
sion on the receding side of NGC 253 at radii of 11.5', 12.8' 
and 19.0', as well as Ha emission on the receding side at a 
radius of 11.5'. The details concerning the identification of 
extended DIG emission are decribed in Section 4. 

No [N 11] emission was seen on the approaching side. The 
3cr upper limit for the EM is ~ 0.09 cm -6 pc (see Section 4). 
The presence of the [N 11] emission on the receding side could 
be the result of a phase jump in the emission line, and the 
detection was just above a 3<r detection. Fig. U shows that 
the filter transmission on the receding side was almost max- 
imal for the [Nil] line, about 80 per cent. However, on the 
approaching side, the transmission of this line drops down 
to 50 per cent, which would make the 3cr detection difficult. 
This could explain why no emission of this line was seen on 
the approaching side. 

For the receding side, no [Nil] emission was seen be- 
tween 12.8' ± 0.7' and 19.0' ± 0.9'. No significant 3<r nor 
even 2a detection was seen. Between these radii, the para- 
sitic emission varied more, i.e. was noisier. Since the [N 11] 
emission that was detected (r = 11.5', 12.8', 19.0') was al- 
ready difficult to detect, being just barely above a 3a level, 
the noisier parasitic emission could explain why no [N 11] 
was detected between 12.8' and 19.0'. Another explanation 
could be that the DIG is in the form of discrete clouds, which 
would mean that emission from it can only be seen at dis- 
crete radii. Even if these clouds are embedded in a diffuse 
background, the EMs from this background might be too 
low to be detected. 

A question remains unanswered. It concerns the Ha 
emission. In the pie-shape procedure, Ha was detected up 
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to a maximum radius of 11.5' ± 0.6', which coincides wit h 
the maximum radius of the H I study of lPuche et all(|l99ll ). 
The strength of the diffuse emission line in this radii bin (re- 
ceding side, r — 11.5') is about a factor of 10 lower than the 
strength of the line in the previous radii bins. No Ha emis- 
sion is observed beyond, on both sides of the galaxy. Fig. 
[4] shows that for Ha, the transmission probability is only 
about 50 per cent for the receding side, and about 80 per 
cent for the approaching side. Hence, detecting Ha emission 
beyond the Hi is much more probable on the approaching 
side with the observing tools of this study, and yet, no 3cr 
detection was seen. 

Bland-Hawthorn et al. (1997) mention the possibility 
that their detections on the receding side could be explained 
by accretion of a faint gas-rich object at large radius. This is 
supported by an image only showing extended stellar light 
at large radii on the receding side (see Bland-Hawthorn et al. 
1997, for more details). This could explain why no extended 
DIG was seen on the approaching side. 

Considering that we were only able to detect emis- 
sion at specific radii and not continuously across all radii, 
the dumpiness of the gas would argue for an accretion 
phenomenon. It has been suggested t hat the extended 
disc in M31 is the re sult of accretion (|lbata et al.l 120051 ; 
IPenarrubia et al.ll2006l ). Also, as mentioned in Section 4.3., 
the [Nil] line we see at 6566.4 A and could be an extremely 
redshifted line with a radial velocity of ~ 850 km s - , which 
would again argue for accretion of a gas rich object onto the 
galaxy. 

However, our observations do not entirely support this 
idea. First, the fact that emission was only seen along the 
major axis argues that this emission could be associated 
with the galaxy. Second, as mentioned earlier, the isophotes 
in the UV GALEX images are symmetric out to r = 13'. 
If there was an accretion phenomenon, the isophotes would 
have been disrupted. This is supported by the isophotes in 
the B band image taken from the Digitized Sky Surveys. 
Images of this survey were taken from the Oschin Schmidt 
Telescope on Polamar Mountain and the UK Schmidt Tele- 
scope. These isophotes also extend until r = 13' and remain 
symmetric in the outer radii. Third, our Ha rotation curve 
is consistent with the Hi rotation curve. Both curves seem 
to behave normally at extended radii, supporting the idea 
that their are no major disruptions. Fourth, the rotation ve- 
locity we obtain for the diffuse Ha emission at r = 11.5 is 
consistent with that of the phase jumped [N n] emission at 
the same radius. The [N n] emission therefore seems to trace 
the same rotation as Ha. If the Ha is not disrupted up until 
11.5', then the [Nil] emission shouldn't be, at least for the 
bin at 11.5'. 

For the reasons mentioned above, we make the assump- 
tion that the extended diffuse gas we detected is truly asso- 
ciated with the galaxy, and that the [N n] emission we see at 
6566.4 A is the result of a phase jump, making its velocity 
(~ 400 km s _1 ) co nsistent with the Ha vel ocity field. 

The study by iKennicutt et all (|200Sft found a mean 
[Nil] /Ha ratio of 0.48 for NGC 253, which includes the 
[N n] lines at 6548 A and 6583 A. However, many stud- 
ies have shown that this ratio can significantly change in 
areas located outside of dense Hn regions and where the 
DIG component is still present (jRevnolds et al.ll 19951 : iRandl 
Il997h . i.e., at radii larger than the neutral component for 



NGC 253. [N n] is a forbidden line, and can therefore be 
seen in areas of diff use gas wher e spo ntaneous emis sion can 
become important. IRandl l| 19971 ) and IRandl (|l998T l showed 
for the galaxy NGC 891 that the [Sll]/Ha, [Nll]/Ha and 
even [Om]/Ha ratio increase significantly with increasing 
distance to the midplane. For the DIG component of the 
Milky Way, the forbidden lines of [N n] and [S n] were also 
found to have intensities with respect to the Ha line of a 
few tenths to unity, and even larger l|Haffner et al.ll2009T l. 
This is much greater than that observed in classical H II re- 
gions. Hence, forbidden lines could be better probes of the 
extended disc. This could explain the detection of [Nil] at 
radii larger than 11.5' rather than Ha. 

The following parag raphs compare the resu l ts fou nd in 
this paper with those of iBland-Hawthorn et all (|l997l ). Ta- 
ble [3] shows the EM values for the diffuse Ha and [N n] 
lines, both for this study and that of Bland-Hawthorn et 
al. (1997). The error of the flux values of our study are esti- 
mated to be 30 per cent. The flux valu es for the diffuse emis- 
sion in IBland-Hawthorn et all (| 19971 ) 's paper are given in 
milli-Rayleigh (mR). The SHASSA survey, used to calibrate 
our data, adopts the conversion that 1R is equal to an EM 
of ~ 2 cm~ 6 pc. We have used this conversion throughout 
this p aper. From this conversion, the lBland-Hawthorn et al] 
( 1997 ) fluxes give the EM values shown in Table[3] The EMs 
fo und in this study are the r efore consistent with those found 
bv lBland-Hawthorn et~aH (119971 ). 

Bland-Hawthorn et al. (1997) used the TAURUS-2 FP 
interferometer at the 3.9m Anglo-Australian Telescope for 
their observations on NGC 253. Four fields, each of 5' each, 
were observed on the receding side of the galaxy. They used 
the FP Staring technique, which gave a single, very deep 
spectrum of a diffuse source for each field of view. The spec- 
trum in a 40 A band is spread out radially from the optical 
axis across the field. This means, that for the spectrum ob- 
tained at each field, an emission line seen at a wavelength 
is an emission line coming from a precise radius in the field. 
They were able to see the Ha and [Nil] (6548 A) line, with 
a total of 10 velocity measurements on their rotation curve. 
Two fields were placed at extended radii - the first was ex- 
posed for 9 hr and led to a [N n] velocity measurement at 
10.5' and a Ha velocity measurement at 12'; the second was 
exposed for 5 hr and led to a [Nil] velocity measurement 
at 13.5' and a Ha velocity measurement at 15'. The veloc- 
ities found for the very extended and diffuse emission are 
shown in Table [3] corrected for the systemic velocity and 
inclinat ion found in this study. 

In iBland -Hawthorn et al] (|l994h . they demonstrated 
that the TAURUS-2 FP interferometer could reach an Ha 
EM of 0.02 cm -6 pc at the 3a level in 6 hr across the full 
FSR (57 A, in their study). Our calculations gave a la up- 
per limit for the EM of ~ 0.03 cm~ 6 pc (see Section 4), 
and we observed the receding side of the galaxy for a to- 
tal of_6Jiours ; _Oirr_s_tudy was not as deep as that 
of IBland-Hawthorn et all (|l994l ). but was deep enough to 
detect extended EMs and went further out. 

According to Table and the bottom-right panel of 
Fig[6 l our velocity measureme nts are consistent with those of 
IBland-Hawthorn et al] l| 19941 ). within the allowed error bars, 
except for two discrepancies. First, no Ha emission was seen 
in our study at 15', whereas IBland-Hawthorn et al] (|l994T ) 
found emission at that radius (EM of 0.08 cm -6 pc). This 
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EM is below our 3er detection and it is possible that our 
detection limits were not sufficient to allow this detection. 

Second, our [Nil] velocity measurement at 12.8' ± 0.7' 
is lying below the rotation curve found by Bland-Hawthorn 
et al. (1997), and more precisely below the [Nil] velocity 
measurement of Bland-Hawthorn et al. (1997) at 13.5'. The 
disagreement between the velocities arises from different ra- 
dial velocity measurements, with a ~ 20 km s~ x difference. 
Although the observations are made with a Fabry-Perot in 
both cases, the technique is not at all the same since there is 
no scanning in Bland-Hawthorn et al. (1997) observations. 
This enables to detect faint fluxes of emission lines with 
much more confidence but is less accurate for measuring the 
velocities since inhomogeneities in the galaxy (such as knots) 
can bias the result because they can show up even if their 
position in the field is not the right one for producing a con- 
structive interference pattern. As a result it seems that the 
error bars on the [N n] velocity points of Bland-Hawthorn et 
al. (1997) are optimistic and underestimated. It is probable 
that the right value of the rotation velocity of NGC253 at 
about 13' radius is intermediate between the value found 
by Bland-Hawthorn et al. (1997) and our value, say around 
210 km/s, hence in continuity with the Ha values found by 
Bland-Hawthorn et al. (1997) at 12' and 15' radius. 

7.2.2 Origin of the very extended DIG 

iHovle fc Ellis] (|l963h mentioned for the first time the exis- 
tence of a diffused matter of free electrons in the vicinity of 
the Galactic plane using a low frequency synchrotron sur- 
vey. Since then, many new properties of the DIG have been 
discovered. Some of these properties inc lude the high ratio s 
of forbidden lines with respect to Ha (|Randlll997l , Il998l ). 
Furthermore, it was found that several edge-on less active 
galaxies show a correlation betw een the DIG in the halo a nd 
the star formation rate (SFR) ()Rossa fc De ttmar 20 03a! bl: 



iMiller fc Veilleuxl 120031 ). The study by iRossa fc Dettmarl 
|2003al ) also suggested that there was a minimum SFR per 
area required to mimic t he disc-halo inter action. However, 
a more recent study from lOev et all (|2007h concluded from 
their SINGS Ha Survey (Survey for Ionization in Neutral 
Gas Galaxies) that the DIG does not seem to correlate nei- 
ther with the Hubble type nor with the total SFR of the 
galaxy. 

Nevertheless, many the ories were proposed to e xplain 
the existence of the DIG. iBland-Hawthorn et al.1 (| 19971 ) 
analysed the case of NGC 253. Several ionization and heat- 
ing sources for the DIG were examined, including the fol- 
lowing: metagalactic ionizing background, compact halo 
sources, compact disc sources (white dwarf population and 
hot horizontal-branch stars), ram pressure heating result- 
ing from shock waves, exotic heat sources such as turbulent 
MHD waves, young stellar disc, dilute photoionization, gas- 
phase depletion from grains, Lyman continuum (Lyc) ion- 
izing photons escaping from OB stars. However, the only 
possible source that could explain the observations was the 
one regarding hot young stars in the central regions of the 
galaxy looking out towards the outer part of the warped disc. 
Ionizing photons could escape from these stars through dust 
scattering, and would be more easily exposed to the outer 
regions because of the warped structure. This explanation is 
in agreement with the detection of both UV emission (prob- 



ing hot stars) out to r = 13' (see Section 5) and diffuse 
ionized emission at r = 11.5', 12.8', 19.0'. 

The calculations bv IBland-Hawthorn et all (|l997h also 
ruled out the source from the metagalactic ionizing back- 
ground. Their study determined that the metagalactic ion- 
izing background was around five times less than that re- 
quired to explain the Ha observations. Their calculations 
also showed that the photoionization model, accounting for 
the Lyc photons escaping from OB stars, could not account 
for the observations, since the observed ratio of [N il] (A 6548 
A) /Ha was slightly higher than unity. This implies that the 
ratio of [Nil] (A 6583 A)/Ha was close to three. The latter 
is twice as high as what was predicted by the photoioniza- 
tion model. The photoionization model also fails to explain 
the lack of He I (A 5876 A) emission. Studies show that the 
ionizing flux necessary to explain the observed [Nil] and 
[S ill of the DIG should allo w the existence of the He I line 
(|Domgorgen fc Mathislll994h . 

Other studies have also shown a correlation be- 
tween the DIG in the haloes of galaxies and their syn- 
chrotron radio-continuum. This might suggest that the DIG 
is correlated with cosmic rays and magn etic fields (see 
iDahlem. Lisenfeld fc Goilal Il995l ). However. ISciamal (| 19951 ) 
discussed that it was difficult for cosmic rays to penetrate 
clouds of gas, and that they should not play an important 
role in generating the DIG. Nevertheless, the Symposium 
paper by Heesen et al. (2009) found a total equipartition 
magnetic field on the order of 7 ^G to 18 /J.G for the disc of 
NGC 253. These values remain high. The magnetic field of 
NGC 253 could therefore play a key role in understanding 
the DIG component. The presence of a disc wind in NGC 
253 was also discussed in the paper by Heesen et al. (2009), 
and could explain the presence of extra-planar H I, Ha, and 
soft X-ray emission in the halo. Finally, it is also important 
to note that the Magellanic streams can not explain the [N n] 
emission found in this study, si nce the radial velocities of the 
two do not correlate spatially (|Mathewson fc Fordlll984h . 

7.3 Consequences of a declining rotation curve 

Declining rotation curves in spiral galaxies are rare. 
However, some have been see n in earlier-type systems 
(ICasertano fc van Gorkomlll99ll). as well as in the edge-on 
Sc galaxy NGC 891 (|Sancisi fc Allenl 1 19791 1. The Gassendi 



Ha survey of SPirals (GHASP) also find a sub-sample 
of galaxies showing rotat i on curves slightly declin ing (see 



iGarrido et alllioO 



„ |2004 [20051 ; lEpinat et al.ll2008l ). More- 
over. lDicaire et al.l (|2008l ) find, for the late- type system NGC 
7793, a decline of ~ 25 per cent « max over half the radius 
range. For NGC 253, IBland-Hawthorn et al.l (119971 ) find a 
decline at the last observed point of the rotation curve of 
~ 10 per cent iWx over the last ~ 10 per cent of r max . This 
study finds a decline of ~ 30 per cent u max over slightly 
more than half the radius range. The only way not to have 
a declining rotation curve would be if the inclination was 
smaller by 30° for the last observed velocity point, which is 
quite unlikely (see Fig. |5J). 

The top panels of Fig. Fig. and Fig. |9] show the 
results for the maximum stellar d isc model, respectively for 
the H i data of lPuche et al.l (|l99ll ). the Ha data without the 
DIG points, and the Ha data with the DIG points. In all 
cases, the observed rotation curve can be reproduced by a 
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stellar disc and a gaseous disc. Hence, no dark halo seems 
to be needed. However, in the central regions of the Ha 
rotation curve, the presence of a bulge can be observed as 
seen by the bump in the curve. To reproduce well this region, 
a bulge component could be added to the maximum stellar 
disc model. 

An advantage of maximum disc models is that a 
unique solution can be determined, whereas many degen- 
eracies exist in multicomponent fits including dark haloe s 
(|Carignan fc Freeman] 1 19851 ; Ivan Albada fc Sancisl Il98fj ). 
Maximum stellar discs also seem to reproduce well the small 
scale irregularities in the rotation c urves, since the stellar 
disc also traces these perturbations (|Freemanlll992i ). How- 
ever, as can be seen in Table 21 the (M/L)+ that best fits 
the maximum-disc model is of ~ 12 Mq/L© for all three 
cases. This value is much too high, especially in the case 
of late-type galaxies such as NGC 253, which causes great 
concern for the va lidity of the mode l. 

The study by ICourteau fc Rbd (|l999h showed that the 
Tully-Fisher relation (TFR) seemed to be independent of 
surface brightness for high surface brightness (HSB) galax- 
ies when plotted asM r versus V2.2 {V2.2 is the velocity for 
an exponential disc at 2.2 scale lenghts). In order to repro- 
duce this feature, their models suggested that discs were sub- 
maximum with Udisc/^totai ~ 0.6. Additionally, their study 
suggested that these sub-maximum discs allowed a smooth 
transition between HSB and low surface brightness galax- 
ies, and that maximum discs for the majority of the HSB 
late-type galaxies are not valid. 

The bottom of Fig. [7] Fig. and Fig. show the re- 
sults for sub - maxim um discs, respectively for the Hi data of 
iPuche et af] |l99ll ). the Ha data without the DIG points, 
and the Ha data with the DIG points. For NGC 253, the 
(M/i)* was fixed to 5.0 Mq/Lq to allow Wdisc/wtotai ~ 0.6. 
For all three cases, the sub-maximum disc seems to agree 
well with the derived rotation curves. The decline observed 
with the DIG points is also well reproduced. The model also 
yields very similar dark halo parameters for the rotation 
curve with and without the DIG component. Hence, the de- 
rived rotation curve for NGC 253 using the DIG points truly 
seems to illustrate a declining rotation curve, and hints that 
the end of the mass distribution has started to be probed. 

Moreover, the sub-maximum disc gives a steeper slope 
in the central areas which seems to reproduce better the 
Ha data at these radii. Table [4] shows that the Ha data 
have a higher central dark halo density than the Hi data, 
which can be explained by the steeper slope of the Ha 
data. However, the dark halo velocity dispersion a seems 
to be fairly the same in all three cases. This parameter is 
a measure of the maximum amplitude of the rotation curve 
jBlais-Ouellette et al.lll999h , and since the maximum veloc- 
ity is essentially the same for the Hi and Ha data, simi- 
lar values are expected. Finally, the (M/L)* in the B band 
expected from stellar population models, using the (B- V) 
given by the RC3 cat alogue (0.85 mag), is ~ 3.46 Mq/Lq 
(Bell fc de Jondl200ll ). This is much more consistent with 
the ratio used for the sub-maximum disc than the one de- 
termined from the maximum disc. 

Our results supports Bland-Hawthorn et al. (1997)'s 
suggestion that the rotation curve is declining in the outer 
parts. While their data points were showing a decline be- 
tween 8 and 9.5 kpc, our outermost velocity point shows 




Radius (kpc) 

Figure 10. Cumulative mass profile (M(< r)) of the dark 
halo. An isothermal sphere was assumed for the dark halo, with 
po = 1.564 Mq pc — 3 and r c = 1.0 kpc. These are the best fit 
parameters of the dark halo for the sub-max disc fit to the data 
that include the diffuse emission. The total mass is in units of 
Mq. 

that the decline seems to continue out to ~ 14 kpc. As seen 
in the mass model analysis, there are two possible inter- 
pretations: 1- either the halo is truncated or 2- the halo is 
highly concentrated (see fig. 7, with r c ~ 1 kpc). Although 
the second interpretation is not favoured in view the mass 
distribution usually seen in late-type spirals were r c is more 
of the order of 5-10 kpc (see e. g. Spano et al. 2008), it seems 
to be the right one since Fig. [TO] (cumulative mass curve of 
the dark halo as a function of radius) shows that the halo 
is not truncated since it is continuously rising. This profile 
was derived using the parameters found with the extended 
Ha rotation curve. 

NGC 253 is the second late-type spiral in the Sculptor 
group with NGC 7793 (Dicaire et al. 2008) that shows a de- 
clining rotation curve in the outer parts. For NGC 253 and 
NGC 7793, the halo maximum radius seem to be ~ 10 kpc. 
This is smaller than the maximum radius of ~ 40 kpc sug- 
gested by a dynamical analysis of the Sculptor Group (Puche 
& Carignan 1991). However, this study was assuming that 
the Sculptor group was virialized which is not necessarily 
the case since the crossing time is of the same order as a 
Hubble time. If the group is not virialized, their estimate 
would then be an upper limit. 



8 CONCLUSION 

Deep FP Ha observations were obtained for one of the Sculp- 
tor Group galaxies, NGC 253. The main conclusions of this 
study are as follows: 

1. The velocity field was obtained using three methods. 
The first lead to an initial velocity field reaching a distance of 
11.5'. The second, in order to extract the global kinematical 
parameters, used a gaussian smoothing function (w\ = 55") 
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and lead to a velocity field reaching 10'. The third used a 
pie-shape smoothing, in order to extract the profiles of the 
DIG component. Typical flux values of ~ 0.1 cm -6 pc were 
reached. 

2. The nuclear region showed the presence of a bar- 
like structure lagging in velocity and showing an increase in 
dispersion. 

3. The observations allowed the detection of the DIG 
component through Ha emission at a radius of 11.5' and 
[Nil] emission at radii of 11.5', 12.8' and 19.0', on the re- 
ceding side of the galaxy. No Ha emission was observed at 
radii larger than the neutral component (11.5'). 

4. From the gaussian smoothed velocity field, the global 
kinematical parameters were derived and yielded a mean 
systemic velocity of 241.2 ± 5.4 km s~ , an inclination i 
of 76.4 ± 3.9° and a position angle PA of 230.1 ± 2.1°. A 
warp is also seen through the variation of the kinematical 
parameters with radius. 

5. A rotation curve was extracted using those global 
kinematical parameters and the initial velocity field. This 
rotation curve is in good a greement with the Hi curve de- 
rived bv lPuche et al.l ( |l99lh . although the slope of the inner 
region is steeper (as expected from the beam smearing of 
the Hi data). The data points resulting from the extended 
diffuse gas (Ha; [Nil] (A 6548 A) were added and showed a 
si gnificant decline although sligh tly different from the results 
of (Bland-Hawthorn et all (jl997r ) . 

6. A mass model analysis was obtained for the Ha rota- 
tion curve with and without the velocities derived for the ex- 
tended DIG emission. The two yielded very similar results, 
both for the maximum stellar disc and the sub-maximum 
stellar model. The declining part of the rotation curve was 
also very well modeled, and seems to truly be declining. 

Deep observations of galaxies can yield very interest- 
ing re sults, as was the case for the study by iDicaire et al.l 
l|2008h which allowed confirmation of a truly declining rota- 
tion curve of NGC 7793. In this study, the emission detected 
in NGC 253 of the DIG confirmed the declining nature of 
the rotation curve. This has profound consequences on the 
study of dark matter properties, and could allow very good 
constraints of the parameters of the dark halo. 
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